---
title: "Monitoreo de Mamíferos y otros vertebrados en la Gran Sabana"
output: 
  flexdashboard::flex_dashboard:
    theme:
      version: 4
      bg: "#F1FBFB"
      fg: "#5B1A18"
      primary: "#D67236"
      navbar-bg: "#FD6467" 
      base_font: 
        google: Prompt
      heading_font:
        google: Sen
      code_font:
        google: 
          # arguments to sass::font_google() 
          family: JetBrains Mono
          local: false
    social: menu
    source: embed
    navbar:
      - { title: "Proyecto en OSF", href: "https://osf.io/6qtep/", align: left }
      - { title: "Repositorio GitHub", href: "https://github.com/jrfep/GS-svocc", align: left }
editor_options: 
  chunk_output_type: console
---

```{r setup, include=FALSE}
library(flexdashboard)
library(dplyr)
library(leaflet)
library(sp)
library(ggplot2)
library(DT)
library(htmltools)
library(units)
#library(plotly)
#library(spocc)
#library(rinat)
#library(mapr)
require(tidyr)
require(raster)
require(chron)
#library(sf)
library(tmap)
#library(lubridate)
#library(calendR)

library(RColorBrewer)
palette(brewer.pal(8, "Set2"))

```

```{r thematic, include = FALSE}
thematic::thematic_rmd(font = "auto")
```

```{r dataread, message=FALSE, warning=FALSE, include=FALSE}

here::i_am("doc/dashboard.Rmd")
GIS.data <- here::here("Rdata","GIS.rda")
if (!file.exists(GIS.data))
  download.file(url="https://figshare.com/ndownloader/files/37547995",destfile=GIS.data)
load(GIS.data)

event.data <- here::here("Rdata","all-events.csv")
if (!file.exists(event.data))
  download.file(url="https://figshare.com/ndownloader/files/42055824",destfile=event.data)
eventos_actualizados <- read.csv2(event.data)
coordinates(eventos_actualizados) <- c("long","lat")
crs(eventos_actualizados) <- grd@proj4string
qry <- over(eventos_actualizados,grd)
eventos_actualizados$grid <- qry$OID_

## differences between both event data frames:
eventos_adicionales <- eventos %>% filter(!species %in% eventos_actualizados$species) %>%
  mutate(fotos=as.character(fotos))


eventos <- eventos_actualizados@data %>% 
  bind_rows(eventos_adicionales)

cam <- camaras %>% 
  group_by(grid) %>% 
  summarise(cam = sum(dias.de.trabajo),
            caz = max(caza.celda)) %>% 
  transmute(grid, caz, cam)

tps_xy <- spTransform(tps, CRS("+proj=utm +zone=19 +datum=WGS84 +units=m +no_defs +type=crs"))
xys <- coordinates(tps_xy)


walk <- tps@data %>% 
  group_by(grid) %>% 
  summarise(walk = n(),
            dates = n_distinct(time))

walk_ms <- data.frame()
for (k in unique(tps_xy@data$grid)) {
  ss <- tps_xy@data$grid == k
  mtz <- as.matrix(dist(xys[ss,]))
  mtz[upper.tri(mtz,diag = TRUE)] <- NA
  walk_ms <- rbind(walk_ms, data.frame(grid=k, dist_walk=sum(apply(mtz[-1,],1, min, na.rm=TRUE))))
}

walk <- walk %>% left_join(walk_ms, by = "grid") %>%
  mutate(region = if_else(grid <151, "Warapata", "Kavanayen"))

gpslog <- spLines(subset(tps,tps@data$track_fid ==0),
        subset(tps,tps@data$track_fid ==1), 
        subset(tps,tps@data$track_fid ==2),
        subset(tps,tps@data$track_fid ==3),
        subset(tps,tps@data$track_fid ==5),
        subset(tps,tps@data$track_fid ==6),
        subset(tps,tps@data$track_fid ==7),
        subset(tps,tps@data$track_fid ==10),
        subset(tps,tps@data$track_fid ==11),
        subset(tps,tps@data$track_fid ==12),
        subset(tps,tps@data$track_fid ==13),
        subset(tps,tps@data$track_fid ==14),
        subset(tps,tps@data$track_fid ==15),
        subset(tps,tps@data$track_fid ==16),
        subset(tps,tps@data$track_fid ==17),
        attr=data.frame(fid=0:14))

event_summary <- eventos %>%
  group_by(grid) %>%
  summarise(
    total_spp=n_distinct(species))

on_camera_event_summary <- eventos %>%
  filter(!camara %in% "RAS") %>% 
  group_by(grid) %>%
  summarise(
    on_camera_spp=n_distinct(species),
    on_camera_events=n())

off_camera_event_summary <- eventos %>%
  filter(camara %in% "RAS") %>% 
  group_by(grid) %>%
  summarise(
    off_camera_spp=n_distinct(species),
    off_camera_events=n())


event_summary <- event_summary %>% 
  left_join(on_camera_event_summary, by = "grid") %>% 
  left_join(off_camera_event_summary, by = "grid")

effort_data <- walk %>%
  full_join(cam, by = "grid") %>%
  full_join(event_summary, by = "grid") %>%
  filter(!is.na(grid) ) %>%
  transmute(grid, 
            walk=coalesce(walk,0L),
            dist_walk=coalesce(dist_walk,0L),
            cam=coalesce(cam,0L),
            total_spp=coalesce(total_spp,0L),
            on_camera_spp=coalesce(on_camera_spp,0L),
            off_camera_spp=coalesce(off_camera_spp,0L),
            on_camera_events=coalesce(on_camera_events,0L),
            off_camera_events=coalesce(off_camera_events,0L)) 


effort_data$bloque <- factor(grd@data$cuadrado[match(effort_data$grid,grd@data$OID_)])

grd@data <- grd@data %>% 
  left_join(event_summary, by=c("OID_"="grid")) %>% 
  transmute(Celda=OID_,
            Bloque=cuadrado,
            #`Distancia recorrida` = dist_walk,
            #`Dias*cámara` = cam,
            `Eventos fotográficos` = coalesce(on_camera_events,0L),
            `Otros avistamientos` = coalesce(off_camera_events,0L),
            `Especies fotografiadas` = coalesce(on_camera_spp,0L),
            `Especies avistadas` = coalesce(off_camera_spp,0L),
            `Total spp. detectadas` = coalesce(total_spp,0L), 
            )

```


Column1{.tabset .tabset-fade data-width=250}
-------
Detecciones con cámara trampa y avístamientos incidentales

### Esfuerzo de muestreo {data-width=245}

```{r esfuerzo}

cuadro_2 <- effort_data %>% 
  mutate(Region = if_else(bloque %in% 1:6,"Warapata","Kavanayen")) %>%
  group_by(Region) %>% 
  summarise(
    `Nr. bloques` = n_distinct(bloque),
    `Nr. celdas` = n_distinct(grid),
    `Nr. celdas muestreadas` = sum(if_else(cam>0 | dist_walk>0,1,0)),
    `Nr. celdas con cámaras` = sum(if_else(cam>0,1,0)),
    `Días*cámara`=sum(cam),
    `Nr. celdas visitadas` = sum(if_else(dist_walk>0,1,0)),
    `Distancia recorrida (km)` = sum(set_units(dist_walk,'m')) %>% set_units('km') %>% round(2),
   .groups = "keep") %>% t()

colnames(cuadro_2) <- cuadro_2[1,]

cuadro_2[-1,] %>% knitr::kable()

```

### Observaciones por especie {data-width=245}
```{r taxalist}


cuadro_1 <- eventos %>% group_by(Especie = species) %>% summarise(`Eventos fotográficos` = sum(!camara %in% "RAS"),`Otros avistamientos`=sum(camara %in% "RAS")) 

  DT::datatable(
    cuadro_1,
    #colnames = ,
    #caption = tbl_caption,
    options = list(
      dom = 't', 
      pageLength = -1,
      order = list(list(1, 'desc'),list(2, 'desc'))),
    #container=my.container,
    rownames = FALSE,
    width = '100%', 
    height = '100%') #%>%
#  formatRound(c("min_RS", "mean_RS", "max_RS", "mean_RS_cor"),digits=3)



```



Column2 {.tabset .tabset-fade}
-------


### Mapa
```{r map, eval = TRUE}
tmap_mode("view")
map1 <- 
  tm_shape(gpslog, name = "Recorrido (GPS)") +
  tm_lines(col = "maroon") +
  tm_shape(subset(grd,Bloque %in% 1:6), name = "Warapata: Unidades de muestreo") + 
  tm_polygons( border.col = "#F1BB7B", alpha=.5) +
  tm_shape(subset(grd,Bloque %in% 7:10), name = "Kavanayen: Unidades de muestreo") + 
  tm_polygons( border.col = "#FD6467", alpha=.5) +
  tm_minimap()

tmap_leaflet(map1)

```


### Caminata {data-width=245}

```{r "walking distance", fig.showtext = TRUE}
ggplot(data=walk, aes(x = walk, y = dist_walk, colour = region)) +
  geom_point() +
  scale_colour_manual(values=c("Kavanayen"="#FD6467","Warapata"="#F1BB7B")) +
  labs(x="GPS log points", y = "Distancia (metros)", title = "Distancia recorrida en cada unidad de muestreo (celda)")
```


### Rastros y otros avistamientos {data-width=245}

```{r xtraKavanayen}
filter(eventos,  camara %in% "RAS") %>% 
  transmute(Especie=species, `Tipo de evidencia`=evento, Bloque=bloque, Celda=grid, Año=ano, Mes=mes,Día=dia) %>% 
  DT::datatable(options = list(
      pageLength = 25,
      order = list(list(4, 'asc'),list(5, 'asc'),list(6, 'asc'))),
    rownames = FALSE)
```


### Esfuerzo de muestreo (expandido)
```{r obsbyyear}

effort_data %>% 
  group_by(bloque) %>% 
  summarise(`nr. de celdas`=n_distinct(grid),
            `Puntos del log del GPS` = sum(walk),
            `Distancia recorrida (km)` = sum(set_units(dist_walk,'m')) %>% set_units('km'),
            `Días cámara` = sum(cam),
            `Eventos fotográficos` = sum(on_camera_events),
            `Otros avistamientos` = sum(off_camera_events)) %>% 
  DT::datatable(options = list(
      pageLength = 25),
    rownames = FALSE) %>%
  formatRound(c("Distancia recorrida (km)"),digits=2)

```

### Información sobre la sesión de R
```{r sessionInfo}
sessionInfo()

```
Column3{data-width=150}
-------

### Warapata

### Eventos fotográficos

```{r}
on_camera <- effort_data %>% filter(bloque %in% 1:6) %>% summarise(sum(on_camera_events)) %>% pull
off_camera <- effort_data %>% filter(bloque %in% 1:6) %>% summarise(sum(off_camera_events)) %>% pull
n_especies <- eventos %>% filter(bloque %in% sprintf("B%02d",1:6)) %>% summarise(n_distinct(species)) %>% pull
valueBox(on_camera, icon = "fa-camera", color = "#F1BB7B")
```

### Otros avistamientos

```{r}
valueBox(off_camera, icon = "fa-binoculars", color = "#F1BB7B")
```
### Especies

```{r}
valueBox(n_especies, icon = "fa-paw", color = "#F1BB7B")
```

### Kavanayen

### Eventos fotográficos

```{r}
on_camera <- effort_data %>% filter(bloque %in% 7:10) %>% summarise(sum(on_camera_events)) %>% pull
off_camera <- effort_data %>% filter(bloque %in% 7:10) %>% summarise(sum(off_camera_events)) %>% pull
n_especies <- eventos %>% filter(bloque %in% sprintf("B%02d",7:10)) %>% summarise(n_distinct(species)) %>% pull
valueBox(on_camera, icon = "fa-camera", color = "#FD6467")
```

### Otros avistamientos

```{r}
valueBox(off_camera, icon = "fa-binoculars", color = "#FD6467")
```

### Especies

```{r}
valueBox(n_especies, icon = "fa-paw", color = "#FD6467")
```